Dynamical exchange interaction between localized spins out of equilibrium 
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The electron mediated exchange interaction between local spins adsorbed on two-dimensional 
surface is studied under non-equilibrium conditions. The effective spin-spin interaction is found to 
depend both on the spin-polarization of the substrate and the excitation spectrum of the local spins. 
For spatially anisotropic spin-polarization of the substrate, the spatial dependence of the interaction 
comprise components decaying as sin(2fcF-R)/(2fcF-R) and sin{2kF R) / {2kF R)'^ . 
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The excitation spectra of spin systems strongly de- 
pends of the type of interactions that are involved. The 
magnetic moment of e.g. single Co [IHS]) Fe [9l ITO] . 
Cr [To], and Mn [11] atoms become strongly anisotropic 
due to symmetry reduction in the interaction with elec- 
tron medium. Studies of inelastic scattering processes 
of layered materials [5] and single atoms [3HZI ISHI3J 
have given deepened insight to the excitation spectra 
of various elements, which then provide further detail 
to the understanding of the involved interactions, phys 
For magnetic systems, the interactions between the lo- 
cal spins can be of different character, which is often 
modeled using e.g. the Ising and Heisenberg Hamilto- 
nians, but also anisotropic models such as e.g. XY- 
or anisotropic Heisenberg Hamiltonians. Regardless 
of model, the interaction parameters describe a phys- 
ical interaction between the spins, which result from 
different mechanisms. The spin-spin interaction may 
be direct in the sense that the exchange Coulomb in- 
tegral {J 'ipl{r)ipl^,{r')V(T,r')'>p„{r')'4>„/{r)drdr') is non- 
negligible, or of indirect nature, e.g. super-exchange 
or double-exchange. Of particular interest is the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction 
[TiHTO] . which is generated by a coupling between the 
local spins and the surrounding electron medium, such 
that the spin-spin exchange interaction is mediated by 
the electronic environment. 

In this paper, we address the electron mediated ex- 
change interaction between localized spin under non- 
equilibrium conditions in two-dimensional systems. The 
question is pertinent to recent measurements using e.g. 
scanning tunneling microscopy (STM) where local non- 
equilibrium conditions are created by the tunneling cur- 
rent. It is demonstrated that the resulting spin-spin ex- 
change interaction depends on the spin-polarization of 
the electron medium and on the excitation spectrum 
of the localized spins. For spatially anisotropic spin- 
polarized surface electrons the spin-spin interaction com- 
prise the Ising, Heisenberg, and Dzyaloshinski-Moriya in- 
teractions, where the latter is shown to asymptotically 
decay as sm{2kpR)/{2kpR). In contrast to previous 
studies of the non-equilibrium RKKY interaction [T7|ll8j. 
we here also include the proper time-dependence of the 



local spins. 

The presence of the local spin (or magnetic) moments, 
results in a spatially inhomogenous surface electron spin- 
polarization which can be transformed into a spatially 
non-uniform spin bias distribution between the spin- 
projections of the surface electrons. Under the spin bi- 
ased conditions, electrons flow between the local spin mo- 
ments, however, different spin projections travel in differ- 
ent directions, thus, establishing a net equilibrium. The 
setup is, hence, reminiscent of the electron spin resonance 
situation discussed in Refs. [TH [20l and is ideal for inves- 
tigating the electron mediated exchange interaction in 
terms of non-equilibrium formalism. 

We begin by considering localized spins Sr at the posi- 
tions r interacting with a continuum, treated in the closed 
time-path Green function formalism [21 which was re- 
cently applied to spin dynamics in a Josephson junction 
[22] and three-dimensional metallic systems [23] . We cal- 
culate the partition function (in units: h = c = 1) 



Z[S,,{t)] =tr Tee 



iS 



(la) 



S =5wzwN + 5cxt + f [Uk + Ht^I, (lb) 
■)dt= / {■)dt+~ / (Odt-, (Ic) 



where we have omitted unimportant contributions from 
the electron gas with quadratic dispersion and isotropic 
effective mass m, as we are considering conduction elec- 
trons in the continuum approximation. The STM tip 
is assumed to have negligible effect on the spin cluster. 

5wzwN = E./S,(i) • [Sr{t) X Sr{t)]dt/Sl, S, = |Sr|, 

is the Wess-Zumino-Witten-Novikov (WZWN) term de- 
scribing the Berry phase accumulated by the local spins. 
The trace runs over the degrees of freedom for the elec- 
trons in the tip and substrate in order to provide an ef- 
fective spin action, which in the present situation repre- 
sents the interaction of the magnetic spins with a non- 
equilibrium environment. S^-Kt represents the coupling 
between the system with the external electromagnetic 
field. The Hamiltonians inside the contour integral de- 
fine the (Kondo) coupling between the local spins and the 
surface electrons, "Hk — —VuJk Ylr Sr(i) ■s(r, t), and the 



coupling to external electrodes T-Lt which generates a tun- 
neling current in and/or out from the two-dimensional 
surface. For example, recent STM measurements mo- 
tivates to model the tunneling current between the tip 
and surface [13] using Ut = Er Epk^a' 4<t('^'^'^'^o + 
Tia„„, ■ Sr)cwe*''-'"+*'^(*) -t-iJ.c, where p (k) denotes the 
momentum for electrons in the tip (substrate), whereas 
To, Ti are the (p and k dependent) rates for the di- 
rect and exchange coupled tunneling. Here, s{r,t) is the 
electron spin density, whereas Vu and Jk defines a unit 
surface element and the Kondo coupling to the electrons. 
(j){t) = e/_^ Vsd{t')dt' gives the energy shift due to the 
bias voltage Vgdit) applied between the tip and surface 
and (T is the vector of Pauli spin matrices. 

The procedure in I22H23J yields the effective action 



S =5„zwN + /"^[5MsB(r,t) +3^^\r,t)/e] ■ Sl{t)dt 
- {vuJk? I Yl Sr(i) • •^'■(r, r'; t, t')Sl,{t')dtdt' 

rr' 

+ ^ / E Sr (i) • J^'' (r, r'; t, t')Sj, {t')dtdt', (2) 

rr' 

where S^{t) = [S{t+) + S(t_)]/2 and S^{t) = 
S(i+) - S(t_), whereas J];,(r,r';t,i') = i~i)0{t - t') x 
([si(r,t),Sj(r',i')]) is the retarded spin GF of the surface 
electrons. j*^^^(r,t) = j*^^)(r,i)z and j^^''(r,r';i, t') are the 
spin-polarized current density and spin current density, 
respectively, between the tip and the substrate generated 
by the spin-imbalance and non-equilibrium conditions in 
the electrodes [Ml US]. This, general, formulation of the 
action is motivated from the perspective of recent tun- 
neling experiments. In this paper the focus, however, is 
on the third term to the right in Eq. ^, which rep- 
resents the RKKY interactions as it emerges from the 
(Kondo) coupling between the localized spin moments 
and the surface electrons. 

Owing to the general non-equilibrium conditions, the 
retarded spin GF is expressed in terms of the lesser and 



< 

greater surface electron GFs G>(r, r'; t, i'), that is, 

X cr^G< (r', r; t\ t) - CT*G<(r, r'; i, t') 
xa^G>(r',r;t',t)], (3) 

where the trace tr 5 is taken over spin space of the sur- 
face electrons. For non-interacting but spin-polarized 
surface electrons we write the real space GFs according to 

G> (R; r) = / G>(k; T)e^'"^dk/(27r)2, where R = r - r' 
and T = t — t' . The lesser and greater forms of the GF 
can be written 

G>(k;r)=^(a° + a.Aa^JG>(k;r)/2, (4) 



where G>(k;r) = (±i)/(±eko-) exp(-iek<TT) with Ckcr = 
£k + (j\A\/2, whereas cr" is the identity matrix, and f{x) 
is the Fermi function. The effective spin-splitting A is 
generated by the surface electrons due to coupling to in- 
ternal and/or external spin degrees of freedom. A can 
partially be due to e.g. the spatially inhomogeneous 
mean field —VuJk X]r(^"") which is generated by the ad- 
sorbed spins, and partially due to e.g. spin-orbit inter- 
actions in the surface, pertinent to recent STM measure- 
ments of Co/Ft (111) [3 . Using Eq. Q and the identity 
(A • (t)(B • ct) == (A • B)crO + i(A X B) • o- [25] we find, 
after some algebra, 

S2(t). J-'-(R;t,i')Sj,(f') = 

= ^ 53F...(R;T){2a^„<,,,[S^(i) • A]K{t') ■ A] 

(T(T' 

+ [l-|A|V^,<,,,]S2(t).SMi') 
-^K,-al,,,]A■[Sl{t)xSl,{t')]]^ (5) 

where the dynamical range functions F„ai are given by 
^,.,(R;r) =(-z)0(r)y'[/(£k'.') - /(£k.)]e^^''-'^')-^ 



X e 



-i(ekCT-ek 



rfk dk' 

(2^(2^ 



(6) 



Below, we shall calculate those functions explicitly. Be- 
fore we proceed, however, we note that the electron medi- 
ated spin-spin interaction described in Eq. ([s]) comprise 
three different kinds of interactions; Ising, Heisenberg, 
and Dzyaloshinski-Moriya (DM) types of interactions, re- 
spectively, in agreement with Ref. [27j . 

The first term (- [S«(i) • A][S^,(t') • A]) is a general- 
ized Ising-type of interaction. That it is of Ising-type can 
be understood since it provides the interaction between 
the spins projected onto the direction of the field A. By 
rotating the reference frame such that e.g. A — Az, 
one finds that [S«(i) • A][S^,(i') • A] = ls.^Sf^''[t)Slf{t'). 
The Ising interaction vanishes for non-spin polarized con- 
duction electrons, since the range functions are spin- 
independent, i.e. Faa' = F, under such conditions. 

The first two contributions are expected to be present 
between spins interacting via metallic or semi- conducting 
medium. The last contribution is, on the other hand, 
expected to arise in anisotropic systems [551 US]- Here, 
this anisotropy is generated by the spin polarized surface 
electrons. For non-chiral spin-polarization of the surface 
electrons, the range functions F^^ = F^^^ which leads to 
that the DM interaction vanishes, as expected. 

It is interesting to note, that the DM interaction is 
non- vanishing whenever the local spins create a spatially 
inhomogeneous spin-polarization of the surface electrons. 
In this sense, the induced spin-polarization can be viewed 
as an effective defect induced spin-orbit interaction. 

Next, we calculate the dynamical range functions. The 
angular integrals in Eq. ^ results in the factors Jo{kR), 



where Jo{x) is the Bessel function. The Fourier transform 
of the exponential e^^'-^'^"^^^'"'''^ is given by (w — Eko- + 
iS)^^{u!' ~ £k'cr' — iS)^^. Thus, for quadratic dispersion 
Ek = fc^/2A'o, iVo = ?7i, changing momentum to energy 
integrations, and carrying out those energy integrals, give 



oo+'i(5 



oo+Z(5 



[U,{u;')-U{co)] 



(7) 

where Hn (x) is the Hankel function, flaa' — |A|((t — 
a')/2, R = RV^No, and /„M = /(w + a|A|/2). Using 
that (-i)0(T)e-^^^ = J{n-x+id)-^e-'^^^dn/2Tr, we can 
finally write F^^,(R;t) = / F^<,-(R, f7)e-'"^df7/(27r), 
where 



i^,,,(R,r!)=H) 



Ni 



[Uioj~n) + Uu;)]H^o'\RV^) 



4 



'2^ 



(8) 



showing that the dynamical range function depends on 
time T, the spin bias ^a-a'j a-nd on the spin-chemical po- 
tential iicr = ep — cr| A|/2. Eq. (Is]) is the central result of 
this paper and below we analyze a few of its consequences 
to the electron mediated spin-spin exchange interaction. 
The static regime, which corresponds to assuming 
frozen spin moments, is defined for f2 = 0. For small 
spin biases Vt^ra' /sf ^ 1, such that H^ ' [R^/uj + ^aa') ~ 
Hq {R^/u))^ the integral in Eq. ([8| can be analytically 

0), for which 



calculated for low temperatures [T 

No 



F.^>(R) 



2n 






Jn{Rka)Yn{RK) 



-\Jl{RK)-Y^{RK 



(9) 



where k^ = \/k-^p — cr| A|7Vo with the Fermi vector kp = 
\/2NoeF, whereas Yn{x) is the Neumann function. Here, 
the real part captures previous results [5DH52] . in the 
spin-degenerate limit k^ — > kp, whereas the imaginary 
part accounts for the retardation and damping effects of 
the interaction due to the electron medium. 

For 2kpR ^ 1, the asymptotic expansion of Eq. ^ 
[50] leads to that ReF^^'(R) - sin{2k„R)/{2k„Rf, in 
agreement with previous results [5UVl32j . which provides 
the usual spatial decay for the isotropic electron medi- 
ated Heisenberg-like exchange, see Eq. ([5]). Analogously, 
the asymptotic expansion of the imaginary part of Eq. 
^ gives ImFo-CT' ~ cos{2k „ R) / {2k aR)^ ■ In contrast, the 
electron mediated DM interaction asymptotically decays 
as sin{2k pR)/ {2k pR), which can be seen from the fol- 
lowing observation. The DM interaction depends on the 
range functions as F^^ — F^^, see Eq. ([S]). By Taylor ex- 
panding the imaginary part of this difference and using 



that fccr ~ fci? — a\A\No/{2kp), one finds that it asymp- 
totically reduces to 



ImVa^,F^;(R) 



|A|iVo sin2A:Fi? 
2kpR ' 



2k p 



(10) 



which to the best of our knowledge has not been re- 
ported previously. Hence, despite the DM interaction 
depends quadratically on the anisotropy field A, it tends 
to dominate over the Heisenberg, and Ising, exchange for 
sufhciently large distances between the spins. Thus, in 
absence of effective magnetic fields acting on the local 
spins, two spins configure themselves perpendicular to 
one another when being separated by a sufhciently large 
distance since the DM interaction dominates their cou- 
pling. A coUinear alignment of the spins is typically fa- 
vorable whenever the spins are close to one another, since 
the Heisenberg interaction provides the strongest contri- 
bution to the coupling. The cross over distance at which 
the DM interaction begins to dominate over the Heisen- 
berg interaction is roughly given by Rq ~ kp/{2\A\NQ). 
This cross over distance is obtained under the assumption 
that No\A\/kp ^ 1, which is reasonable from the point 
of view of Ref. [33], where kp ~ 0.17 A~^ whereas the 
spin-splitting due to the Rashba effect -\/]Vo|A| ~ 0.012 
A~^ (|A| ~ 1 meV). In terms of these values, the cross 
over distance Re ~ 1, 180 A. For larger spin-splitting, 
the above approximations are not valid, however, we ex- 
pect that the effect of the DM interaction will be even 
stronger and that the cross over distance is significantly 
shorter. 

For the remainder of this paper, we consider the lo- 
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FIG. 1: Energy (fi) and spatial {R) dependence of the real, 
(a), (c), and imaginary, (b), (d), parts of _F-f|(R, f2). The 
dotted lines in panels (a) and (b) indicate the energies at 
which the line scans in panels (c) and (d) are extracted. In 
panels (c), (d), the plots are vertically shifted for clarity. Here, 
kp = 0.17 k-\ fitJ. = 10 meV at T = lOK. 



cal spin moments to be time-dependent, i.e. f2 7^ 
0, as they would be under the influence of e.g. an 
effective magnetic field. In Eq. (|5|, the dynam- 
ical range function F^^i (R, r) mediates the interac- 
tion between the spins signified by S'^^(i') and S^^)(i), 
which have different time- arguments. Using that e.g. 
S(i)(t) ^ JS^^\x)e'''*dx/{2'!T) and integrating over t 
and i' [c.f. Eq. (pi)], we can write, for instance, 
the Heisenberg type of exchange interaction between 
two spins located at r and r' as —{vuJK)'^J2cra'i^ ^ 

analogously for the other two types of interactions. The 
exchange interaction between the localized spins, hence, 
strongly depends on the excitation spectra of the spins. 
In Fig. [1] we plot the real, panels (a), (c), and the imag- 
inary, panels (b), (d), parts of i^cro-' (R, fi) for a given 
spin bias, Q^^ ~ 10 meV, and temperature, T = 10 K. 
As expected from previous results, and from the above 
discussion, in the regime near i7 = 0, F^a-' acquires an 
oscillatory decaying behavior as function of R. 

The expression given in Eq. (|8| suggests to inter- 
pret the effective exchange interaction as the interfer- 
ence between (spin-dependent) charge density waves with 
their frequencies set by the spin-chemical potential ^a, 
the spin bias flaa', and the excitation spectra of the 
localized spins. Particularly, the waves are described 
in terms of the wave vectors k = \/27Vqw and k^^r' — 
y2]Vo(ar^TT+TIov), respectively. Fig. |l](a), (b), dis- 
play how the oscillatory character of F^a' (R-, ^\ as func- 
tion of i?, changes as the energy fl varies along the verti- 
cal axis. The period of the oscillations is shorter for en- 
ergies ri ^ compared to the period in the static regime. 
In the static regime, the frequency of the charge density 
wave differ only by the spin bias ^aa' , which typically is 
small compared to the Fermi energy. Hence, the charge 
density waves are almost in phase with one another. For 
finite energies, fi ^ 0, the difference in frequency be- 
tween the two waves increases, such that the waves go 
out of phase. Hence, due to the incommensurability of 
the waves, the resulting period of Fa-a' (R, ^), as function 
of R, changes for increasing \il\. 

In conclusion, we have studied the electron mediated 
spin-spin exchange interaction under non-equilibrium 
conditions for localized spins embedded in a two- 
dimensional system. It was demonstrated that the range 
function depends dynamically on time, the spin excita- 
tion spectrum, and the spin-bias between the spin chan- 
nels in the electron medium. This leads to that the elec- 
tron mediated exchange interaction between the local- 
ized spins is determined by the spin-polarization of the 
electron medium as well as of the excitation spectra of 
the local spins. In the case of spatially anisotropic spin- 
polarized surface electrons, the electron mediated spin- 
spin exchange interaction comprise the Ising, Heisenberg, 
and Dzyaloshinski-Moriya type of interactions, capturing 



the static case previously reported p7j . It was, moreover, 
shown that earlier results for the range function, which 
were derived for systems in the static regime [57[ I5D1 - 
l32| . can be straightforwardly extended to slowly fluctu- 
ating spins. Particularly, the Dzyaloshinski-Moriya in- 
teraction was shown to decay as sm{2kpR)/{2kpR) for 
weakly spin-polarized electrons. 

The author thanks A. V. Balatsky, A. Bergman, O. 
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for support from the Swedish Research Council. 
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